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We investigate the formation of fuUerene-like structures from hot Carbon gas using classical molec- 
ular dynamics, employing Brenner's potential. In particular we examine the influence of different 
annealing strategies on fullerene yield, which is characterized by the distribution of coordination 
numbers and polygon numbers. It will be shown that the fullerene yield strongly depends on the 
annealing strategy. Furthermore, we observe a close relation between polygon formation and the 
number of atoms surrounded by three atoms. 
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I. INTRODUCTION 

A. Fullerenes and Nanotubes 

Since their discover jffl in 1985, fullerenes have been the 
subject of intense studies, both experimentally and the- 
oretically. These clusters have been identified as roughly 
spherical carbon "cages" . The fullerene family includes 
molecules ranging from the (probably unstable) C20 to 
fullerenes which contain hundreds or even thousands of 
carbon atoms. The most stable fullerenes are the soccer 
ball-like Ceo (also known as buckyhall) and the C70. 

All fullerenes have in common the same principal 
structure: They are composed of sp^-hybridized carbon 
atoms, therefore being closely related to graphite. How- 
ever, since they are not planar, they can't be made up 
entirely of hexagons. The spherical structure is typically 
achieved by insertion of pentagons. Also Carbon nan- 
otubed^, discovered 1991, have a similar structure ~ the 
main difference to fullerenes is the enormous aspect ratio 
(length vs. diameter). 

Nanotubes can be interpreted as roUed-up graphite 
sheets, where the precise procedure of roUing-up decides 
about the electronic structure of the nanotubes (metallic 
or semiconducting). So in the absence of defects, nan- 
otubes are composed entirely of hexagons - except the 
caps which close the tubes. 

In the language of classical geometry, fullerenes can be 
described as polyhedrons, so Euler theorem applies: For 
an arbitrary polyhedron with V vertices, E edges and 
F faces, the Euler characteristic C = V — E + F is an 
invariant, for all simple polyhedrons, C = 2. From this 
theorem it can be deduced that any structure which is 
topologically equivalent to a sphere and only composed of 
riQ hexagons, ns pentagons and nj heptagons must fulfill 

715 - = 12 (1) 

with arbitrary ng. So all fullerene molecules will contain 
exactly 12 pentagons (unless there are heptagons present 
to compensate for a higher number of pentagons). There- 
fore, C20 is indeed the smallest possible fullerene, com- 
posed of 12 pentagons and no hexagons at all. 



Since pentagons induce stress, structures which con- 
tain adjacent pentagons show reduced stability [isolated 
pentagon rule, IPR). Obviously, for less then 60 atoms 
there is no possibility to avoid this occurrence, and also 
for 61 < n < 69 atoms no structure exists which fulfills 
the IPR. In perfect icosahedral Cgo two pentagons are 
always separated from each other by one hexagon, this is 
the most stable configuration, followed by the rugby-ball 
like C70 which again fulfills the pentagon rule. 

Further increasing the number of atoms does not en- 
hance stability. Multiple adjacent hexagons show a 
graphite-like structure that is more severely disturbed 
by the pentagons the bigger the molecule gets. So the 
yield of those fullerenes (for example at arc discharges) 
is again much smaller than that of Cgo or C70. 

B. Molecular Dynamics Studies 

A natural tool to study the dynamics of fullerene sys- 
tems is molecular dynamics (MD), where the equations 
of motion for the nuclei are integrated numericall}'^ The 
most important difference between various types of MD 
simulations is the way to implement the interaction be- 
tween the nuclei. 

Depending on whether t he ele ctronic structure is simu- 
lated by ab-initio method^^^^, tight-binding methodd^ 
or by more or less elaborate force fields, the accessible 
system sizes and simulation times vary dramatically. 

While there is a wide variety of force- field approaches, 
certain forms are particularly well-suited for the simu- 
lation of (hydro-)carbon syste ms: Bu ilding on the ap- 
proach of Abell and Tersoff" iHBEU^ Donald W. Bren- 
ner has developed a potential^Sn^ for hydrocarbons that 
is closely related to the potentials obtained in the em- 
bedded atom methods (EAM). It contains empirical pa- 
rameters which were fitted to reproduce the properties of 
hydrocarbons as well as those of carbon in diamond and 
graphite. 

Its first use was the study of chemical vapor deposi- 
tion, but found widespread use in various simulations of 
carbon-based materials. This potential has been incor- 
porated in the BrennerMD cod^^^, which was used for 
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the simulations of sec. HT] 

MD studies of fuUerenes by now encompass a broad 
range of topics, among others 

• studies of nanotube growtbPIMIIIlIIlSl^ 

• charge structure of nanotube growth from small 
cluster^, 

• nanotube cap formation and energetic^^ll 

• reactive collision of graphite patches^^ which lead 
to formation of the nuclei of nanotube, nanocage 
and other structures. 



• collisions between fullerene^^MU and between 
fuUerenes and other molecule^^H^ 



• the impact of fuUerenes on graphite materials^, 

• formation of diamond crystallites inside nested car- 
bon fuUerenes.''^ 

A more comprehensive overview has been given else- 
wher^. Most interesting in the context of the present 
article are those studies which deal with formation and 
fragmentation of fuUerenes or fuUerene-like structures. 

FuUerene stability and fragmentation mechanism have 
been investigated with a variety of computational meth- 
odd^, in particular with tight-binding methods. There 
Ceo turned out to be stable against spontaneous disinte- 
gration for temperatures up to 5000 KP^l. Similar results 
have been obtained for other Cn-systems! ^'^ * '^^ ! ^ . Most 
common fragmentation products are dimers, rings and 
chains; the fragmentation temperature first increases lin- 
early with cluster size, but becomes nearly constant for 
fuUerenes composed of 60 or more atoms. 

Other simulations'^'' find structural changes in Cgo and 
C70 between 3000 and 4000 K and bond breaking around 
5000 K. For the melting and evaporation of C20, C'go and 
C240 one finds"^* a transition from the low-temperature 
solid to a floppy "Pretzel phase", melting at about T = 
4000K and conversion to carbon chain fragments when 
increasing the temperature to T « lO^K. 

A notable feature of disintegration is the ejection of 
small cluster^^H (in particular C2, sometimes also C3) al- 
ready below the melting temperature. Similar results for 
Desintegration and cage formation have been obtained 
using the TersofP^ or Brennei'^ potential. 

When, studying fullerene formation and using mono- 
cyclic, bicyclic, and tricyclic rings as precursors'^, the 
most efficient temperature for the formation of the cage- 
shaped structure from a ring is about 3000 K. Also the 
growth process of higher fuUerenes through adduction 
of small carbon clusters and carbon atoms has been 
studied^ In this case small clusters and single atoms 
easily adsorb on the surface of fullerene cages which have 
defects, on collisions at thermal velocities; during an an- 
nealing process, the attached clusters are soon incorpo- 
rated into the network of the fullerene cages. 



Summarizing these results, most studies find that dur- 
ing the annealing process, one has a sequence of typical 
steps: 

• For very high temperatures, one has a gas of single 
atoms or dimers. 

• When cooling down, the dimers tend to form chains 
that stay short first, but grow to length of 0(10) 
atoms below T = 7000K. This phase ("Pretzel 
phase"'^ can be interpreted as liquid-like.^ 

• Below T Ki 5500K most chains begin to form 
vaguely spherical structures where most atoms 
establish an additional bond to reach carbon's 
favored coordination number of three. In partic- 
ular, hexagons, which are sign of locally graphite 
respectively fuUerene-like structures, form. 

• The precise transition temperature for the change 
to fuUerenes of course depends on the model 
employed for interaction, but the growth and 
removal of defects seem to be most efficient around 
T ^ 3000 K. 

The resulting structures are typically graphite- or 
fuUerene-like with lots of defects. The number of de- 
fects can be reduced by slow annealing, and it has been 
possible (e.g. by isolating a fullerene precursor) to obtain 
perfect fullerene molecules in some simulations. 

All these simulations still have a severe drawback, the 
times accessible by these simulations are still short com- 
pared to the timescales relevant for most physical pro- 
cesses relevant for the formation of fuUerenes, which is 
supposed to happen at 0(s). Simulations, however, can 
reach in the case of rather small systems with severely 
restricted types of force fields at most 0(100 ns). In 
the present paper we address the question, whether this 
shortcoming can be compensated by the choice of an ap- 
propriate annealing strategy. 



C. Organization of the Paper 

The paper is organized as follows: After the introduc- 
tion to fuUerenes (sec. lAl and corresponding molecular 
dynamics simulations (sec. I B I , we turn to the basic sim- 
ulation setup, described in section |II] 

Our analysis tools are described in section |III| 
where we discuss distribution of coordination number in 
sec. |III A| and polygon numbers in sec. |IIIB"| In sec. |IV| we 
present and discuss our results; in sec. |V| we summarize 
them and give a brief outlook. 



^ Note that this phase is not always founds. We suspect that this 
is an artifact of the force field implementation, as for example a 
fixed neighbour list. Such a constraint makes the formation of 
long chains extremely unlikely since (in the case of 60 atoms) for 
a given dimer only four out of 58 atoms are capable of attaching 
to this dimer. 
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FIG. 1: Annealing hot carbon gas: snapshots. This simulation is similar to the ones performed with annealing strategy Ti, 
explained in sec. |ll]and illustrated in fig. |3] One can clearly recognize both fuUerene-like and chain structures. 



II. SIMULATION SETUP 

Even with modern computing power, the times ac- 
cessible to MD simulations are still short compared to 
the timescales relevant for most physical processes to be 
studied by this method. This also applies to the forma- 
tion of fuUerenes, which is supposed to happen at 0{s) 
while simulations are typically limited to 0(100 ps), since 
one elementary MD step corresponds to 0{fs). At best 
(studying small systems with cheap and thus severely re- 
stricted types of force fields) one can reach 0(100 ns). 

One might expect that this shortcoming can be par- 
tially compensated by the choice of an appropriate an- 
neahng strategy, e.g. if repeated heating-coohng cycles 
can significantly improve the fuUerene yield - and this is 
the main question to be studied in this article. 

A. Examples for Annealing Simulations 

We have simulated the annealing of carbon gas from 
T « lO^'K to room temperature. The simulations have 
been performed using the BrennerMD cod^^, which 
makes use of Langevin dynamics, employing the Fungi- 
mol graphical user interface^. 

An example for such type of simulation is shown in 
fig. [l] (4024 carbon atoms in a box of 100 A x 100 A x 100 A 
with a simulation time of t — 245ps). The cooling was 
done in decrements of 100 K with 5000 thermalization 
timesteps of 0.5 fs length in between. So the total sim- 
ulation lasted ttotai = 245 ps plus a few femtoseconds of 
initialization steps without thermostat. 

The typical stages discussed in sec. |IB| are clearly vis- 
ible here: the dimer stage, the liquid-like chain phase, 
the emergence of triple-bound atoms which tend to form 
planar structures and thus the formation of fuUerene-like 
structures. 

Of course many defects occur even in the graphite- or 
fuUerene- like structures; a significant amount of atoms 
is still contained in long chains even after cooling to 
room temperature. These chains have (at least for 
even numbers of atoms) higher ground-state energy than 
ring-based structures (including graphite, graphene and 
fuUerenes), but they seem to form easier and are favor- 
able at high temperatures due to entropic reasons. 

During the annealing process they should convert to 
smaller rings, but here this is often impossible because of 
the rapid cooling process. There are some similarities to 
the formation of an amorphous state via rapid cooling. 
If - in a very similar setup - the monotonic decrease of 
temperature is replaced by repeated heating an cooling 
cycles, the final state typically looks like the one illus- 
trated in fig. [2j 

One still has a considerable amount of long chains, but 
now they are wound up and enclosed by roughly fuUerene- 
like graphitic cages. The resulting structures can be in- 
terpreted as partially frozen droplets with a solid shell, 
but a liquid core. 
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FIG. 2: Cut through two final structures obtained by repeated 
heating and cooling (for two simulations starting with a gas 
of A'' = 1620 atoms) 

B. Setup for the Main Simulation 

In order to characterize the infiuence of annealing 
strategies in a systematic way, we have repeatedly per- 
formed several simulations, employing different annealing 
strategies. In all these simulations the system consists 
720 carbon atoms, initially distributed on a simple cubic 
lattice in a box of 80 A x 80 A x 80 A. 

All simulations start with some initialization steps at 
^Initial = 9900 K and end at the same final tempera- 
ture Tfinai- In all cases, the total simulation time was 
t — 300 ps (600 000 timesteps of 0.5 fs). Data were col- 
lected at the beginning of the simulation and every 5 ps 
thereafter. Seven different annealing strategies have been 
employed, which are also depicted in fig. |3] 

(Ti) After an initialization period of 24 000 steps (corre- 
sponding to a simulation time of 12 ps) at TinitiaiK, 
the temperature is lowered by AT — 100 K every 
6 000 timesteps (corresponding to 3ps) until Tfinai 
is reached. 

(T2) After the same initialization period as in T2, the 
system is rapidly cooled {AT — —100 K every 1 000 
steps) down to T = 6000 K at which the system re- 
mains for 97000 steps. This process (rapid cooling 
by lOOOK in 5ps, followed 48 ps at constant tem- 
perature) is repeated until T = 2000 K is reached; 
afterwards rapid cooling is continued down to Tfinai- 
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(Ta) After the same initialization period as in T2, the 
system is rapidly cooled (AT = —100 K every 1 000 
steps) down to T = 3000 K. After that, the temper- 
ature is rapidly (AT — ±100 K in 0.5 ps) raised to 
T = 6000 K and lowered again down to T 3000 K. 
This heating-cooling cycle is repeated in total eight 
times; afterwards the system is further cooled down 
(AT = -lOOK in 0.5ps) to Tfinai- 

(T4) After an initialization period of 56 000 steps (28 ps), 
the system is rapidly (AT = — lOOK in 0.5 ps) 
cooled down to T = 3 000 K, heated to T = 5 800 K, 
cooled to T = 2 800K, heated to T = 5 600K, 
cooled to T = 2 600 K and so on, until, after heat- 
ing from T = 1 600K to T 4400K the system is 
cooled down to Tfinai- 

(T5) After an initialization period of 23 000 steps, the 
system is rapidly cooled (AT = —100 K every 1 000 
steps) down to T = 4000 K. After that it is cooled 
slowly (AT = — 50K every 7 000 steps) down to 

Tfinai- 

The last two strategies are just variations of (T3), but the 
heating-cooling cycles operate between T = 2000 K and 
T = 5000 K for (T6) respectively between T = 1000 K 
and T — 4000 K for (T7) . The annealing curves for theses 
seven strategies are displayed in fig. [3] 



III. ANALYSIS TOOLS 

In this article we focus on two characteristic quantities: 
the distribution of coordination numbers and types of 
polygons present in the system. 



A. Coordination Numbers 

The coordination number z (number of direct neigh- 
bors) is characteristic for the phase of carbon. A gas 
of carbon atoms has z — Q, for carbon dimers one has 
z = \. Carbon chains and isolated rings have z = 2. For 
the sp^-hybridized carbon atoms in graphite, graphene, 
nanotubes and fuUerenes one has z = 3, while diamond 
structures are characterized hy z — A. 

The bond length of sp^-hybridized carbon in graphite 
layers is 1.42 A, for sp^-hybridized atoms in diamond one 
has I.54A. The carbon-carbon bond lengths in Cgo are 
I.40A for adjacent hexagons and 1.45A in pentagons. 
Similar values one has to expect in other fuUerenes or 
nanotubes. Therefore setting up a threshold of a* := 
I.6OA for the carbon-carbon bond seems to be a reason- 
able choice. In our simulations, bonds are attached to 
all pairs of atoms that have a distance smaller than the 
threshold a*. While other values for this threshold have 
some influence on the results, this effect is small within 
reasonable limits. We will use a normalized coordination 
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FIG. 3: Annealing strategies Ti to Tj as described in sec.|II B| 



number count 



N ' 



where is the number of atoms with coordination num- 
ber z and N is the total number of atoms. 



B. Number of Polygons 

Also the number of polygons present in the system 
is characteristic. FuUerenes are mostly composed of 
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hexagons and several pentagons, whereas pure carbon 
nanotubes as well as graphite and carbon fibers are made 
from hexagons only. Heptagons may break the convexity 
of fuUerene structures, they repeatedly appear as defects 
during cooling processes. 

Simulations will be characterized by the normalized 
number of polygons 



number of polygons of type a 
iV/2 



The polygon type can be pentagons (a = p), hexagons 
(a = h), or heptagons {a — 7). The number of polygons 
is normalized to the graphenc limit of hexagons, which is 
=y if there are N atoms. 

In order to determine the number polygons of length 
L, we use a ring number algorithm which is built on a 
tree-search as depicted in fig. |4] 

The idea is as follows. For each atom the algo- 
rithm checks whether it is part of a ring of length L 
(L = 5, 6, 7). The algorithm is repeated for the various 
values of L and for all atoms. Starting from the initial 
atom, the list of nearest neighbor (n.n.) atoms is deter- 
mined. The n.n. list contains those atoms which are not 
further away from the reference atom than the threshold 
for the bond length. 

In turn, for each atom in the n.n. list the corresponding 
n.n. list is generated. The procedure is repeated L times. 
The resulting structure can be regarded as a graph with 
atoms as nodes and bonds as edges. 

So, roughly speaking, for each particle all possible 
paths of length L starting from the selected atom and 
leading in each step to one of the nearest neighbors of 
the current atom, are determined. From all these paths, 
those for which starting and end point coincide (closed 
paths) are selected. 

Next we eliminating those structures, in which bonds 
are repeatedly traversed (as shown in fig. [5^ and|5]3). 
When all particles have been checked, the resulting num- 
ber of detected closed paths of length L equals 2L times 
the number of polygons in the system, since each atom 
of a polygon is used as initial atom and for each atom 
there are two possible paths, one clockwise, the other 
one counter-clockwise. 



IV. RESULTS AND DISCUSSION 
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The goal of this paper is to reveal the details of the 
polygon formation in order to develop an improved an- 
nealing strategy. To this end, MD simulations for the 
different annealing strategies presented above are per- 
formed and analyzed. In all simulations, the total num- 
ber of atoms, initial and final temperatures, as well as 
the total number of MD steps is the same. 

In figs. [6|i and [6]3 the final results are summarized in 
form of the normalized coordination number count and 
the normalized number of polygons for the various an- 
nealing strategies. An interesting observation of fig. 



FIG. 4: Tree-based polygon search: basic algorithm; graphical 
representation of six-, seven- and eight-atomic paths 

is the fact that the ratio of the number of hexagons to 
pentagons and heptagons is almost the same in all cases. 
An exception is strategy T3, in which the polygon yield 
is negligibly small in any case. 

A comparison of fig.|6^ and|6]D reveals that the polygon 
yield is roughly proportional to the number of atoms with 



7 




(a) 
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FIG. 6: Results: (a) normalized coordination number count 
for z = (0,1,2,3,4) after annealing runs Ti to Ty; (b) 
number of polygons (pentagons p, hexagons h, heptagons 7), 
normalized to the graphene limit, for annealing runs T\ to T7 



coordination number 3. This implies that polygons are 
not formed as isolated rings, but rather immediately as 
locally connected planar structures. The ranking of the 
annealing strategies in view of the polygon yield is T^, 
Ti, T2, T5, T4, T7, Tq. The most ineffective strategy of 
all is clearly T3 with almost no polygons, but long chains 
instead. 



The expression 

number pentagons — number heptagons 

^Mi = ^ (2) 

could be expected to characterize the number of 
fuUerenes present in a system, provided all carbon atoms 
are precisely bounded to three other atoms. 

In our simulations, however, we find the number of 
heptagons to be typically larger than the number of pen- 
tagons. At second thought, this should not really be sur- 
prising since in the "chain phase" the formation of longer 
rings is more likely than the formation of shorter ones. It 
shows, however, that we are still quite far from a system 
which only contains fuUerenes. 

According to fig.[6]the strategy T3 has a very poor poly- 
gon yield, while the opposite is true for the best strategy 
Tg, although the two annealing strategies look fairly sim- 
ilar. The only but crucial difference is that in the first 
case the heating-coohng cycles operate between 3000 K 
and 6000 K, while in second case the temperatures are 
2000 K and 5000 K. There are two effects which are ex- 
pected to be responsible for that: 

• In the temperature range of (2000 K, 3000 K) lies 
the crossover temperature T3 , at which the number 
of atoms with coordination number 3 start to dom- 
inate over those with coordination number 2. In 
this range the formation of graphite- or fuUerene- 
like structures seems to be most efficient. 

• At the same time the heating-cooling cycles of Tg 
do not exceed (contrary to the T3) the onset tem- 
perature To of polygon formation. As we see in 
fig. |9] the system reacts almost instantaneously on 
temperature changes, therefore, heating above Tg 
would break up the elementary building blocks of 
the polygons. 

Fig. [7] contains the the detailed time evolution of the 
normalized coordination count and the normalized num- 
ber of polygons during the anealing process for strategy 
T2, in which a stepwise cooling is performed with five 
intermediate temperatures, such that the system is al- 
most in thermodynamic equilibration at the intermediate 
temperatures. Based on the results in fig [7^ we define re- 
gions in which atoms of a specific coordination number 
dominate. We conclude that atoms with coordination 
number 2 dominate already at T = 10000 K down to 
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FIG. 7: Annealing strategy T2: (a) temperature progression; 
(b) normalized coordination number count; (c) normalized 
number of polygons. 



FIG. 8: Annealing strategy T5: (a) temperature progression; 
(b) normalized coordination number count, (c) normalized 
number of polygons. 



T — 6000 K. From T — 6000 K down to a crossover tem- 
perature T3* = 2000 K, atoms with coordination number 
2 predominate, while below T3* atoms of coordinations 
number 3 are in the majority. In the observed tempera- 
ture range, we found almost no atoms with coordination 
number 4, the prerequisite for the diamond structure. 

The onset of polygon formation, according to fig |7]3, 
takes place at about tp = 90 ps, while the forma- 
tion of atoms with coordination number 3 is roughly at 
to = 125 ps. The polygon formation happens right af- 
ter the cooling from T = 6000 K to T = 5000 K, while 
the formation of atoms with coordination number 3 takes 
place later at the same temperature. Hence, the polygons 
are formed first as ringlike structures before they combine 
locally to form planar structures. Interestingly, the op- 
posite is the case if the system is cooled more rapidly 
without reaching thermodynamic equilibrium, as is ob- 
served in fig. |8] 

The onset of polygon formation in fig [8] is at a signifi- 
cant lower temperature than in fig. [7] due to the inertia 
of the system. Interestingly, the simulation time before 
polygon formation sets in is the same in both cases. The 
comparison of coordination number count in both strate- 



gies reveals that the system passes through comparable 
phases, i.e. the depletion of isolated atom, the forma- 
tion of dimers, the increasing and decreasing number of 
atoms in chains, followed by the gradual formation of 
atoms with coordination number 3. 

The best of all scrutinized strategies is Tg, where the 
systems passes through several heating-cooling cycles. 
The results are depicted in fig. [9j Although the fraction 
of atoms with coordination number 3 seems converged in 
fig. [7]i for T 2000^ at a value of q,=3 w 0.5, the an- 
nealing in Tq leads to a significant higher yield of about 
0.7 when T = 2000-?^ is reached during the last cool- 
ing cycle. At the same time the number of atoms with 
coordination number 2 is reduced by the same amount. 
Similarly, the number of polygons is roughly increased by 
40% in strategy Tg, when T = 2000 K is reached during 
the last cooling cycle. 

Exemplary simulations with a strategy similar to Tq, 
but a larger number of heating-cooling cycles yield the 
results given in table [ij The extension from eight to 16 
cycles significantly increases the number of triple-bonded 
atoms (mostly at the expense of double-bonded ones) and 
the number of hexagons. A further extension to 24 cycles 



9 



(a) 'pj'j^j Temperature progression -type 6: 

lOOOOy-;- 

8000- \^ 
6000- \, 

4000- \^ >/\ /'x. /x y'\ 

2000- \/ \y \X \ 

0.0 50.0 100.0 150.0 200.0 250.0 300.0 

t[ps] 

(b) 



0.0 100.0 20O.O 3CO.0 

normalized coordination number count: ""[P^ 

no neighbours 
— — one neighbour (diniers) 
w*Mw two neighbours (chains) 

lAAA three neighbours (graphite/graphene/fullerenes) 
^esB four neighbours (diamond) 



(c) 




M-M-M pentagons, hexagons, i>o-D heptagons 



FIG. 9: Annealing strategy Te: (a) temperature progression; 
(b) normalized coordination number count, (c) normalized 
number of polygons. 





8 cycles 


16 cycles 


24 cycles 
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0.013 


0.008 


0.010 


92 


0.243 


0.178 


0.171 


93 


0.731 


0.789 


0.796 


94 


0.013 


0.025 


0.024 


9p 


0.119 


0.083 


0.078 


9h 


0.358 


0.444 


0.442 


9r 


0.152 


0.122 


0.128 



TABLE I: Effects of longer simulation runs on the fuUerene 
yield for strategy Ta- The simulations have been performed 
for a specific intial configuration and eight, 16, 24 instead of 
eight) 



only marginally changes the results, so for this particular 
strategy we seem to have reached a plateau. 



V. SUMMARY AND OUTLOOK 



We have simulated the structures resulting from an- 
nealing a hot carbon gas with different strategies, i.e. 
different forms of temperature evolution. The struc- 
tures have been analyzed with respect to the coordination 
number count and the number of pentagons, hexagons 
and heptagons. 

We have found the yield of fullerene-like structures to 
be quite sensitive to the strategy employed. Thus the 
fuUerene yield can either be significantly enhanced or al- 
most reduced to zero by appropriate annealing strategies. 

The best strategy included several heating-cooling cy- 
cles, operating below the fuUerene disintegration temper- 
ature (corresponding to the onset temperature of polygon 
formation), but encompassing the crossover temperature 
between double- and triple-bonded atoms 



Since timescales in MD simulations are short compared 
to those of most physical processes, an optimized an- 
nealing strategy can be an important means to obtain 
more realistic results for MD simulations of fuUerene for- 
mation. Alternatively, repeated heating-cooling cycles 
which, performed in the optimal temperature window, 
gave the best fuUerene yield, might also occur (on larger 
timescales) in nature and enhance the fuUerene yield e.g. 
in soot. 

While some scrutinized strategies were obviously supe- 
rior to uniform annealing, it is not at all clear whether 
we have found an even close-to-optimal strategy in or- 
der to maximize the fuUerene yield. Here a stochastic 
optimization (employing simulated annealing or genetic 
algorithms) of the annealing strategy might provide valu- 
able information on how far the limits can be pushed for 
short-time simulations. 
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